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Using a momentum representation of a magnetic von Neumann lattice, we study a two- 
dimensional electron in a uniform magnetic field and obtain one-particle spectra of various periodic 
short-range potential problems in the lowest Landau level. We find that the energy spectra satisfy 
a duality relation between a period of the potential and a magnetic length. The energy spectra con- 
sist of the Hofstadter-type bands and flat bands. We also study the connection between a periodic 
short-range potential problem and a tight-binding model. 
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A physical system sometimes shows a peculiar fractal structure of the spectrum if the system has two scales of 
periods. In the nearest-neighbor (NN) tight- binding model with -a magnetic flux, Hofstadter obtained one-particle 
spectra of butterfly-shape and discovered a multi-fractal structure.u Experimentally, it may be a challenging theme to 
observe the spectruni-jO|£|the fractal structure. Actually, interesting phenomena have been observed in experiments of 



lateral super lattices.H □ The system with the cosine potential has been studied well theoretically.El In the case of tha 
antidot array, however, the potential term includes a large number of the cosine functions with different Jivave lengths 
Tbe steeper the potential of the antidot is, the harder to solite the eigenvalue problem is. SilberbauerQ and Kiihn et 
alu studied a steep antidots potential problem. Huang et alu solved the periodic short-range potential problem in a 
finite system numerically. Analytic properties of the problem are unknown. 

In the present paper, we study eigenvalue problems of the two-dimensional electron systems defined on an infinite 
plane with the periodic short-range potentials in the lowest Landau level (LLL) analytically and numerically. We find 
Hofstadter-type bands and flat bands. The former bands satisfy a duality relation. Namely, the spectrum at one value 
of t is connected with that of where t is the flux $ penetrating the unit cell of the periodic potential normalized 
by the piiit flux $o, i.e. t = $/$o- Our duality is concerned with flux and is different from the Aubry and Andre's 
dualityjl3 which is concerned with hopping strengths in the NN tight-binding model. We also study the connection 
between a periodic short-range potentiajUjjroblem and a tight-binding model. 

The magnetic von Neumann latticald'ta is a representation of the quantum system in a uniform magnetic field and 
has quite useful properties in studying the system with periodic potentials. In particular, its lattice structure varies 
and we can select the-suitable one in accordance with the potential lattice. Using a momentum representation of the 
von Neumann lattice,LJ we give a proof of the duality and show energy spectra in the periodic potential problems. 

In a two-dimensional system under a uniform magnetic field -B, the von Neumann lattice basis is formed by the 
direct product of harmonic oscillator eigenfunctions |//) of relative coordinates (^,77) and coherent states \amn) of 
guiding center coordinates {X,Y). The coherent states are defined by 

[X -f iY)\a,nn) = a2m«|amn), (1) 

where a = 2^1% / eB , Zmn — mujx + ncjy {m,n £ Z) and uix,ujy are complex numbers that satisfy Im[a;*Ci'j,] ~ 1.0 
For our purpose, it is convenient to use the Fourier transformed basis as 

where /3(p) is a normalization factor defined by 

/3(p) = (2Imr)^ e^T,Pli)^(Pl±IP^\r), (3) 

and T = —u)x/u!y For t = i, the von Neumann lattice becomes a square lattice. For r = e'^'^/'^, it becomes a 
triangular lattice. The Fourier transformed basis is an orthonormal set of extended states and obeys a nontrivial 
boundary condition 

|;,p-f27rN)=e-'*(f'^)|/,p), (4) 

where (p{p, N) = 7r(iV^ + Ny) - Nyp^. 

If the magnetic von Neumann lattice has a periodicity commensurate with the periodicity of the external potential 
V'(x), the one-body potential problem becomes easy to treat. This happens when t is equal to q/p with relatively 
prime integers p, q. Let us consider an arbitrary regular lattice of short-range potentials of t = q/p: 
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y(x) = a^Va^5'''^\z + aN^quj,^ + a^ujy). (5) 

N ^ 



The matrix element of the potential in the LLL is given by 



(0,p|F(x)|0,p') (6) 

X (27r)2<5(p^ - + 27r(^ + N,))5{p'y - Py + 2nNy)e'1'^P' 

where Vp = 0, 1, — 1 and Vq = 0, 1, g — 1. Here, we study the eigenvalue equations in the LLL. The eigenvalue 
equation becomes 

D'^DtP = etp, (7) 

where a p x q matrix D is defined by 

{D{p))r^r, = (}*{p, - 27r(^ + (8) 

p q 

and D'fD is a g X g Hermitian matrix. The magnetic Brillouin zone (MBZ) is the region where \px \ < ir/q and \py \ < tt. 
Each band has a p- fold degeneracy. Consequently, the fundamental region of the MBZ is the region where \p,j;\ < T^/pq 
and \py \ < TT. 

We should note that the rank of D''D is min(p, q) generally. If p < g, the band splits into p subbands and one flat 
band which corresponds to the zero modes of D. The number of the zero modes is q — p. By a linear transformation, 
■0' — Dip, the eigenvalue equation becomes 

DD^' = et/j'. (9) 

DD^ is a pxp Hermitian matrix. Equation (0) is equivalent to Eq. (^) except for the zero modes. Apparently, Eq. (H) 
is obtained by taking the complex conjugate and interchanging p and q of Eq. (^. Therefore, there exists a duality 
relation between two problems of t = q/p > 1 and t = p/q < 1. Actually, the energy spectra E obey the relation 

m - jEil), (10) 

except for the flat bands. If p > q, Eq. (^) is also obtained. Thus, the duality relation is proved in an arbitrary 
regular lattice of short-range potentials. The self-dual point is < = 1 and the physical quantity has critical behavior 
near the point. 

Before solving the eigenvalue problems numerically, we clarify the connection between a periodic short-range 
potential problem and a tight-binding model in the LLL. Let us suppose a square lattice potential V{x.) = 
Vbexp(27ri(ma; -I- ny)/b), where t — {b/af' = q/p. The potential energy term in the second quantized form becomes 



where 



Cr,(p) = &o(p. +7r- 2^^r„py)e-*"f'-' (12) 
q 

and blip) is the annihilation operator of the state |^, p). Equation ( pT| ) is equivalent to the tight-binding Hamiltonian 
with the following hopping term: 

Vhop = Voe^^^^™^~'"^^^'^^"^~"^^^^^''^^™^^'"^^^"^^"^\ (13) 

where m2 — mi = m and n2 — ni — n. This hopping term has a flux per unit cell of a square lattice 27r/t and the 
hopping strength varies with t. If {m,n) = (±1,0), (0, ±1), Eq. (O) becomes the NN tight-binding Hamiltonian. In 
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a square lattice of short-range potentials, the summation is taken over all (m, n) with an equal weight. Therefore, the 
result becomes a tight-binding Hamiltonian with finite-range hopping terms 

i ^ c^(m2, n2)Vhop{m2,n2; mi, ni)c(mi, ni). (14) 

The hopping range is about \/t and the number of relevant terms increases linearly with t. 

In a regular lattice, the one-particle spectrum is given by the solution of Eq. (Q). Since the matrix element of D 
is given in the analytic form, it is easy to solve numerically. Figure 1 shows the spectra for the square and triangular 
lattices. The points at = correspond to the original LLL. There are two marked structures in these figures : 
Hofstadter-type bands and a large gap above the flat bands in t > 1. The origin of the large gap is easily understood 
in a dilute potential limit t — > oo. In this limit, the potential approaches one short-range potential. Its spectrum 
consists of a bound state trapped at the potential and a flat band. The bound state has the energy E — Vq and the 
flat band has the energy E = 0. In a dense potential limit t ^ 0, the potential approaches constant Vo/t. Therefore, 
the asymptotic form of the spectrum in t ^ and t ^ oo except for the zero modes is given by 

E{t)^Vo{l + j), (15) 



which satisfies the duality relation (10) 



To check the duality, we replace E with Et for t < I and replace t with 2 — l/t otherwise in Fig. 1. Then the 
figures become symmetric with respect to t = 1 owing to the duality. The results for the square lattice are shown 
in Fig. 2. The duality is clear in this figure. The patterns in Figs. 1 and 2 are very similar to that of the NN and 
next-nearest-neighbor (NNN) tight-binding model on the square and the triangular lattice. However, in contrast with 
the NN or NNN tight-binding model, a periodicity with respect to the flux does not exist in our model. Instead, the 
duality between t and l/t does exist. 

As an example of an irregular lattice, we study a honeycomb lattice potential numerically. A honeycomb lattice 
can be regarded as a sum of two triangular lattices. This property causes the following eigenvalue equation at 
2t = (flux per hexagon)/<I>o — q/p- 



{DlDi+DlD2)^ = e^, (16) 



where Di and D2 are defined by 



Di^D{p,,Py), D2=D{p^- —,py + ^). (17) 

bmce both d\Di and dJDz are positive definite, the zero modes must satisfy both Diip = and = 0. If 

q — 2p > 0, the zero modes exist and the LLL splits into 2p subbands and one flat band. Otherwise, the LLL splits 
into q subbands. This system does not satisfy the duality relation (|l^). However, the asymptotic behavior of the 
one-particle energy spectrum is also given by Eq. ( [l5| ) from arguments similar to that of the regular lattice. The 
numerical result of the spectrum is shown in Fig. 3. The spectrum also has a Hofstadter-type structure and a large 
gap above the flat bands. 

Next we study the possibility of observing the band structure obtained above. The gap above the flat bands may 
be observed in the magnetoresistance experiment .□ We assume that the antidot potential has a height Vq and an area 
of the base Tq. Vq in Eq. (||) is related to Vq and as a^Vb = TqVq. The magnitude of the large gap above the 
flat bands is of the order of Vq and the correction of the finite size effect of the antidot to the energy is estimated 
as VoO(ro/6) where 6 is a lattice constant of the antidot array. Therefore the gap above the flat bands should be 
observed in ro ^ h. In the current experiments,u r^/h ~ 0.25, which is enough to observe the gap. Taking the Landau 
level mixing into consideration, the eigenvalue equation for the potential of Eq. (^ has the form 

y[Ei5u' + ^I?Wt^(0]^(0 ^ ,^(0. (18) 

Because of the (q— p)-zero modes of i?^'^ for p < q, the l-th Landau band has a flat band with the energy Ei. Therefore, 
the flat bands remain flat in the presence of the Landau level mixing effect in the short-range periodic potential 
problem. The correction of the Landau level mixing to the energy above the large gap is estimated by the second- 
order perturbative calculation as Sll ~ tf^)^/^'-^c- From Vq ^ -Ell, we obtain the condition hujc ^ Vo(?'o/a)^ for the 
gap to survive. Using the realistic valuesB'Q m = O.G7me and Vq — 0.3[meV], this condition becomes 27r?i^/mr§ ^ Vq- 
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This is satisfied in the current experiments. So far, conditions for observation are satisfied. However, the parameter t 
in current experiments, e.g. B = 5[T], h = 200[nm] and t k, 50, is too large to observe the large gap, the Hofstadter- 
type bands, and the critical behavior near t = 1. Thus we hope that the experiment will be made in a finer lattice of 
the antidot array. 

In summary, we solved various periodic potential problems in the LLL using a momentum representation of the von 
Neumann lattice. In a periodic array of short-range potentials, the energy spectrum has three remarkable structures, 
(i) The spectrum has a Hofstadter-type structure, which is commonly seen in the periodic potential problems in a 
magnetic field, (ii) The duality relation exists for a regular lattice potential of Eq. (^). A honeycomb lattice is 
an irregular lattice and obeys a duality relation asymptotically, (iii) There is a large gap above flat bandalj that 
comes from zero modes of the matrix D of Eq. (||). These structures are universal in the periodic short-range potential 
problem in the LLL. The conditions for the observation of these structures were obtained. In addition, the equivalence 
between two problems of a square lattice of short-range potentials and a tight-binding model with a inverse flux was 
shown. 
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(b) TRIANGULAR LATTICE 
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Figurel: Energy spectra of the periodic short-range potential for (a) the square lattices and (b) the triangular lattices. The 
horizontal variable t is the flux per unit cell in the flux unit. Vb is the strength of the potential. 
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Figure2: Energy spectra of Fig. 1(a) are deformed to show the duahty clearly. The vertical variable is tE/Vo for t < 1 and 
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E/Vo otherwise. The horizontal variable is 2 — 1/i for i > 1 and t otherwise. As a result, the spectrum becomes symmetric 
with respect to t = 1. 
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FigureS: Energy spectra for the honeycomb lattices. The horizontal variable t is the flux per half of the hexagon in the flux 
unit. Vb is the strength of the potential. 
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